A note on Wiener-Hopf factorization for Markov 
Additive processes 

Przemyslaw Klusik * Zbigniew Palmowski''' 
June 5, 2009 

Abstract 

We prove the Wiener-Hopf factorization for Markov Additive pro- 
cesses. We derive also Spitzer-Rogozin theorem for this class of pro- 
cesses which serves for obtaining Kendall's formula and Fristedt repre- 
sentation of the cumulant matrix of the ladder epoch process. Finally, 
we also obtain the so-called ballot theorem. 

1 Introduction 

The classical Wiener-Hopf factorization of a probability measure was given 
by Spitzer (1964) and Feller (1970), and has a strong connection to ran- 
dom walks. This result was generalized by Rogozin (1966), Fristedt (1974) 
and other authors using approximation based on discrete time skeletons. 
Greenwood and Pitman (1980) use direct approach which relies on excur- 
sion theory for reflected process. For details see Bertoin (1996) and Kypri- 
anou (2006). Presman (1969) and Arjas and Speed (1973) generalize Spitzer 
identity into different direction, to the class of Markov Additive processes in 
discrete time (see also Asmussen (2003) and Prabhu (1998)). Later, Kaspi 
(1982) proves Wiener-Hopf factorization for a continuous time parameter 
Markov Additive process, where Markovian component has a finite state 
space and is ergodic. The main weakness of the fluctuation identity given by 
Kaspi (1982) is that they involve distribution of the inverse local time, which 
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is seldom explicitly known. Dieker and Mandjes (2006) investigate discrete- 
time Markov Additive processes and use an embedding to relate these to a 
continuous-time setting (see also Breuer (2008) and Rogers (1994)). 

This paper presents Wiener-Hopf factorization for a special, but none 
the less quite general, class of Markov Additive Processes (MAP). For this 
class of processes we give short proof of Wiener-Hopf factorization based 
on Markov property and additivity. We also express the terms of Wiener- 
Hopf factorization directly in terms of the basic data of the process. Finally, 
we derive Spitzer-Rogozin theorem for this class of processes which serves 
for obtaining Kendall's formula and Fristedt representation of the cumulant 
matrix of the ladder epoch process. We also present the ballot theorem. 

The paper is organized as follows. The Section [2] introduces basic defi- 
nitions, facts and properties related with MAPs. In Section [3] we give the 
main results of this paper. Finally in Section U] we prove all theorems. 



2 Preliminaries 



2.1 Markov Additive Processes. 

Before presenting main results we shall simply begin by defining the class of 
processes we intend to work with and their properties. Following Asmussen 
and Kella (2000) we consider a process X{t), where X{t) = X'^^^t)+X'^^\t), 
and the independent processes X^^\t) and X^'^\t) are specified by the char- 
acteristics: qij,Gi,ai,ai,h'i{dx) which we shall now define. Let J{t) be a 
right-continuous, ergodic, finite state space continuous time Markov chain, 
with Z = {!,... ,N}, and with the intensity matrix Q = {qij). We denote 

the jumps of the process J{t) by {Tj} (with Tq = 0). Let {Un^} be i.i.d. 
random variables with distribution function Gi{-). Define the jump process 
by 

n>l iel 

For each j G I, let X^ (t) be a Levy process with the Levy-Khinchine expo- 
nent: 



logE(exp{iaX-'(l)}) = ^j{a) 

= -iaja-\ h / {1 - e'"'^ + ia\y\l\y\<i)vj{dy), 

oo 



oo 



where /^(l A \y\'^)i'j{dy) < oo. By X^'^\t) we denote the process which 
behaves in law like X^{t), when J{t) = i. We shall assume that the afore 
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mentioned class of MAPs is defined on a probability space with probabili- 
ties {Pj : i G I}, where Pi(-) = P(-| J(0) = i), and right-continuous natural 
filtration F = {J^^ : t > 0}. In fact we can consider more general MAP 

where additional jumps uji^ appearing during the change of the state of 
J{t) could also depend on the state J(r„_|_i) (so called anticipative MAP). 
This could be done by considering the vector state space and the modified 
governing Markov process J on it. If each of the measures i^i are supported 
on (—00, 0) as well as the distributions of each U^^^ then we say that X is 
a spectrally negative MAP. These definition and more concerning the basic 
characterization of MAPs can be found in Chapter XI of Asmussen (2003). 

2.2 Time reversal 

Predominant in the forthcoming analysis will be the use of the bivariatc 
process { J,X), representing the process {J,X) time reversed from a fixed 
moment in the future when J(0) has the stationary distribution tt. For 
definitiveness, we mean 

J{s) = J{{t - s)') and X{s)=X (t) - X {{t - s)-) , 0<s<t 

under = X^jgx^*^*- Note that X is also Markov Additive process. The 
characteristics of (J, X) will be indicated by using a hat over the existing 
notation for the characteristics of (J, X). Instead of talking about the pro- 
cess (J, X) we shall also talk about the process (J, X) under probabilities 
{Pj : i el}. Note also for future use, following classical time reversed path 
analysis, for y > and s <t, 

¥i{G{t)eds,-I{t)edy\J{t) =j) 

= Pj {G{t) e ds, S (t) - X{t) e dy\J{t) = i) , (1) 



where I{t) = info<5<t X(s), S{t) 
X{s) = S{s)}, G{t) = sup{s < t 
explain the last identity). 



= supo<s<tX(s) and G{t) = sup{s < t : 
: X{s) = I{s)}. (A diagram may help to 



2.3 Ladder height process 

We start from recalling the representation of the local time given in Kaspi 
(1982) in formula (3.21). For MAP we say that state i € {l,...,iV} is 
regular when ¥i{R = 0) = 1 for i? = inf{i > : t e M}, where M is a. 
closure of M. = {t > : X{t) = S{t)}. Denote by {Un} the stopping times 
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at which R{t-) = and R{t) > for the R{t) = inf{s > t : S{t) = X{t)} - 1 
and J(t) is irregular. Denote 

g ^(Un on X{Un) = SiUn), 

" \ OO otherwise. 

By Theorem 3.28 of Kaspi (1982) (see also Maissoneuve (1975)), for the 
MAP we can define the ladder height process: 

{iL-Ht),H{t) = X{L-\t)),J{L-\t))),t > 0} 

choosing the local time: 

L{t) = L\t)+Y,\{J{Un))ef\ (2) 

S'„<t 

where L'^{t) is a continuous additive process that increases only on M. and 

In) 

e:\ are independent exponential random variables with intensity 1, 

A(i) = E, [(1 - e-^] . 

Obviously to make functional ([2|) measurable we enlarge probability space 
to include the exponential random variables. One can easily verify that 
{L~^{t),H{t),J{L-^{t))) is again (bivariate) MAP (see Kaspi (1982, p. 
185)). For each moment of time we can define the excursion: 

^ r X{L-\t-) + s) - X{L-\t-)) for L-\t-) < L-^t) 

\ d otherwise, 

where 5 is a cementary state. Let Ci^t) = L~^{t) — L^^{t—) be the length of 
the excursion. From ([2]) it follows that the excursion process {{t, et), t > 0} 
is (possibly stopped at the first excursion with infinite length) marked Cox 
point process with the intensity n{J{L~^{t—)),de) depending on the state 
process J{L~^{t—)). Denote by £ the u-field on the excursion state space. 

2.4 Spectrally negative Markov Additive process 

Letting QoG(a) = {qijGi{a)), where Gi{a) = E (exp(aC/^*'')) , for spectrally 
negative MAP we can define cumulant generating matrix (cgm) of MAP 

x{ty. 

F(a) = QoG(a) +diag(V'i(a),...,V'7v(a)), a E M+ . (3) 
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Perron-Probenius theory identifies F (q) as having a real-valued eigenvalue 
with maximal absolute value which we shall label k (a) . The corresponding 
left and right 1 x A'" eigenvectors we label v (a) and h(a), respectively. In 
this text we shall always write vectors in their horizontal form and use the 
usual to mean transpose. Since v (a) and h (a) are given up to multiplying 
constants, we are free to normalize them such that 

V {a) h {aj^ = 1 and vrh {a)^ = 1 . 

Note also that h (0) = e, the 1 x vector consisting of a row of ones. We 
shall write hi (a) for the i-th element of h (a) . The eigenvalue k (a) is a 
convex function (this can also be easily verified) such that k (0) = and 
k' (0) is the asymptotic drift of X in the sense that for each z G X we have 
lim^ioo E(X(t)| J(0) = i,X(0) = x)/t = k' (0). For the right inverse of k we 
shall write ^. 

It can be checked that under the following Girsanov change of measure 



(ml 



._ g7AW-M7)t for 7 such that k (7) < oo (4) 



the process (X, Pj) is again a spectrally negative MAP whose intensity ma- 
trix F^ (a) is well defined and finite for a > —7. Generally for all quan- 
tities calculated for P'^ we will add subscript 7. Further, if F^ (a) has 
largest eigenvalue k,^ (a) and associated right eigenvector h-y (a), the triple 
(F^ (a) , (a) , (a)) is related to the original triple (F (a) , k (a) , h (a)) 
via 

F^ {a) = Ah (7)""'" F (a + 7) Ah (7) I and k,^ (a) = K{a + 7)-^ (7) , 

(5) 

where I is the N x N identity matrix and 

Ah (7) := diag {hi (7) , h^ (7)) • 

Similarly, the time reversed process X{t) is the spectrally negative MAP 
with the characteristics F, h, k. To relate them to the original ones, recall 
that the intensity matrix of J must satisfy 

Q = A^iQ^A^ , 

where A^ is the diagonal matrix whose entries are given by the vector vr. 
Hence according to ([3]) we find that: 

F(q) = A;iF(af A^. 
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Moreover, k {a) = k (a) and A^h (a)^ = v (a)^ (see Kyprianou and Pal- 
mowski (2008) for details). 

The spectrally negative MAP is easier to analyze since its ladder height 
process {L^^{t), H{t), J{L~^ (t))) has explicit matrix cumulant generating 
matrix S(q', a). Let 

r+ := mi{t > : X{t) > a}, 

where a > 0. Denote the generator of the Markov process {J{t^),(i > 0} 
by A(g) on P*(5)_ it solves equation: 

F$(,)(-A(g)) = 0; (6) 

see Pistorius (2005) and Ivanovs et al. (2008). Note that ladder height 
process can be identified as {{t^ , X{t^) = a,J{T^)), a > 0}. It is a 
bivariate Markov additive process with the cumulant generating matrix: 

S{q,a) = Ah($((7))(^('?)I-A(g))Ah($(g))-^ + a 

= Ah($(g))(($(g) + «)I-A(g))Ah($(g))-^ (7) 

for a, q> 0. Above could be deduced from the equalities: 

and the Theorem 1 of Kyprianou and Palmowski (2008) stating that 

E(e-«^^"; T+ < e,; J(r+)) = E {e-^"^^^^^ l^^^+^^y, J(r+)) 

= Ah(^(g + 0)e-(*(^+^)'-'^(''+^))" Ah(«5(g + 0)"' 

= exp {-Ah($(g + + 01 " Mq + e))Ah($(g + 0)"'^} -(9) 

3 Main results 

As much as possible, from now on, we shall prefer to work with matrix no- 
tation. For a random variable Y and (random) time r, we shall understand 
E(y; J(t)) to be the matrix with (i,j)-th elements Ej(y; J(t) = j). For 
an event, A, P{A] J(t)) will be understood in a similar sense. Here and 
throTighout we work with the definition that is random variable which 
is exponentially distributed with mean 1/q and independent of (J,X). Let 
Iij{q) = Fifl{J{eq) = j), in other words 

I{q)=qiql-Q)-'. 
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From now on we assume that none of the processes X* are downward sub- 
ordinators and compound Poisson process. To include compound Poisson 
process X^'^\t) in Theorem[Ui) on the event {J(G(eg)) = i} it is necessary 
to work with the new definition G_{t) = inf{s < t : X{s) = I{s)} instead the 
previous one. 

Theorem 1. (i) For a general MAP the random vectors {S{eq),G{eq)) 
and {S{eq) — X{eq),eq — G{eg)) are conditionally on J{G{eq)) independent. 
Hence for a G R, ^ > 0, 



E 



E e*"^^'^')-^^^^'); J(G(eg)) 



(10) 



A^^E 



JaI{eq)-(Gieq). 



and 



E 



E 



JaI{eq)-i(G{eq). 



'W(G(e,)) 



A-^E 



JaS{eq)-i^G{eq) . 



^V(G(e,)) 



(11) 



A^ 



(ii) For the spectrally negative MAP and a,(, >0, 



E 



-aSieq)-(G{eq). 



;J(G(e,)) =H(g + e,a)-Miag(S(g,0)I(g)e^), (12) 



E 



-a5(e,)-CG(e,) 



H(g + e,a)-'H(g,0)I((?), 



(13) 



E 



'e'^n^,)-iG(^,).j^G{eq)) 



(14) 



q {{q + e)I - F(a))-i A'^Siq + -a) ^ diag ( Siq, 0)I(<7)e 



-1 



E 



,a7(e,)-eG(e,). 



J{eq 



q {{q + 01 - F(a))-i A'^iq + ^ -a)^ [s(g, 0)-^ 



(15) 



Remark 2. Applying Theorem [^i) to the reversed process derives similar 
conclusion for the infimum functional. Namely, processes {{X{t), J{t)),0 < 
t < G{eq)} and {(^(^(e,) +t)- X{G{eq)), J{G{eq) + t)),t > 0} are condi- 
tionally on J{G{eq)) independent. 



7 



Remark 3. For N = 1 (hence A{q) = 0, l{q) = 1) above theorem gives 
well-known identities for the spectrally negative Levy process: 



E 



E 



-a5(e,)-{G(e,) 



^(9 + + «' 



where il){0) = —'^{—iO) is a Laplace exponent of X. Finally, for ^ = 
above theorem gives already known identity for spectrally negative MAP 
(see Kyprianou and Palmowski (2008)): 

E(e°^(-'');J(e,))^(F(a)-9l)T 

= gAv(a>(g))[a($(g)I-A(Q))-i-I]Av(<l>((z))-\ (16) 

which was derived using Asmussen-Kella martingale. 

We prove also the following counterpart of Spitzer-Rogozin version of 
Wiener-Hopf factorization and the Fristedt theorem: 

Theorem 4. Assume that the matrix E ex.p{iaX (1)} has distinct eigenval- 
ues and that for any t,s > 0; 



E 



JaX{t) 



= E 



i{x(t)>o} ;<-'"(*) 



E 



l{X(s)<0}; J{s) 



JaX{s) 



l{X(s)<0}; <^(s) 



E 



JaXit) 



J{t). (17) 



Then 



E 



-a5(e,)-eG(e,). 



exp 



dt 



[0,oo) 



-^t—ax 



and 



E 



exp 



d.t 

^(-00,0) 



-^t+ax 



- l) t-^e-«*P(X(t) G dx; J(t))| l{q) 

- l) t-ie-«*P(X(i) G dx; J(i))| l{q) . 
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The assumption ()17p is satisfied for example for Markov modulated 
Brownian motion X{t) = a{J{t))B(t), where cr is a positive function. 

The following generalization of the Kendall's identity and the ballot the- 
orem also hold. 

Theorem 5. Consider the spectrally negative Markov Additive process X{t). 
If there exist distinct qi, q2, ■ ■ ■ , qn such that vectors h(<I>(gi)), h($(g2)), 
. . . , h($(gAr)) are independent, then 

tP(r+ G dt; J{t))dx = xV{X{t) £ dx; J(t))dt . 

Theorem 6. Let 

X{t) =ct- a{t), 

where {a{t),t > 0} is a Markov Addidive Subordinator without drift com- 
ponent. Under the assumptions of the Theorem O the following identity 
holds: 

F{X{t) e dx,I{t) = 0; J{t)) = ^P{X{t) G dx-J{t)) . 

One can straightforward check that the assumptions of the Theorem [5] 
are satisfied e.g. for X{t) = ct — J{t)N{t), where N{t) is a Poisson process 
and J{t) is a two-state birth-death process. 

In total theorems given here might be seen as a fundamental of the 
fluctuation theory for the (spectrally negative) MAP and might serve for 
the deriving counterparts of well-known identities for the Levy processes. 



4 Proofs 

4.1 Proof of Theorem [1] 

(i) Sampling MAP process {X{t),J(t)) up to exponential random time e^ 
corresponds to the sampling the marked Cox point process (double Pois- 
son point process) of the excursions up to time L(eq). Moreover, since 
conditioning on realization of the process J{t) the point process (t, ej) is a 
non-homogeneous marked Poisson process, we know that, conditioning on 
J(L-i((T^-)) for 

= inf{t >0:eteA}, 

the point process {{t,et),t < a^} is independent of e^A. 
Consider now 

r fL-Ht) ] 

ai = mi it > : l{x(s)=5(s)} ds> 
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and 



a2 = inf{t > : C(et) > e*}, 
where e* is the independent exponential random variable with intensity q if 



1 

et ^ d and e* = 5 otherwise. Note that a2 is for the ^ = {C(e) > e^}. 



If £72 < 0"i, then conditioning on J(L "^((fJi A (T2) — )) = J{L ^(o"2 — )) the 
process 

{(t, et), t < cJi A CT2 and / 5} (18) 



is independent of e, 



0"2 



^cri Acr2 • 



If fJi < (72, then e. 



5 and is also 



independent of the process (fTSj) . Hence conditioning on J(L ^(cri A 0"2 — )) 
the excursion e(jiA(T2 is independent of the process (fTSll . Note also that 

(19) 



G(e,) = L-i((ai A as)-), 5(e,) = H{{ai A as)-) 



and the last excursion e(jjACT2 occupies the final eq — G{eq) units of time in the 
interval [0, e^] and reaches the depth X{eq) — S{eq). This completes the proof 
of the first part of the Theorem [T][^i) . Note that (e^ — G{eq^X{eq) — S{eq)) 
has the same law like (G(eg), / (e^)). The second part of of the Theorem [TJi) 
follows now from the first part applied for the reversed process. 

(ii) To prove Theorem [T||i) we follow Bertoin (2000). Fix n G N setting 



where 



i+ = [nS{eq)]/n, (20) 
stands for integer parts. Applying the strong Markov property at 



time and using ^ yields: 



E 



■,J{eq 



00 



k/n 



-ak/n^ 



; k/n < S{eq) < {k + l)/n; J{eq) 



k=0 
00 



A;=0 



n(I - e-=('?+^'°)V«] 



-1 



nP 



(21) 



(22) 



Taking n — > 00 we have —i- S{eq) and r — >■ G{eq) and hence the left hand 
side of above equation converges by the dominated convergence theorem. 
Thus also right hand side converges. Note that for any matrix A 

1 



-A/n 



n 



-A + o(l/n). 



(23) 
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Thus 



E 



-5G(e,)-aS(e,). 



J(eg) =H(g + e,a)-^B 



for some matrix B using fact that matrix {^{q + C) + ck)I ~ ^{q + is 
invertible for q > 0. Taking = a = we obtain 

B = S{q,0)I{q) 

which completes the proof of (jl3p . Similarly, from (j2ip . 



e »n ;J(r+)=j 



fc=0 



E 



e 



P <r+ ; J(e,) 



and therefore 



E 



-?G(e,)-aS{e,). 



;'/(eq 



E 



-^G(e,)-a5(e,). j^j^^g^^y 



C , (24) 



for 



-1 



P(eg<r+„;J(eg)) 



^ = „l™odiag(^n^eg<r+„ 

= hm P(J(eg)|5(eg) < 1/n) . 

n—^oo 

Now (|12p will follow straightforward from (j25p and ()13p . since by 



(25) 



P (e, < r+„; J(e,)) = 1(g) - P (e, > r+„; J(e,)) = (l - e-^^"'")!/") 1(g) 



and hence 



C= [diag(H(g,0)IeT)] 'H(g,0)I(g). 



(26) 

The identity ()14p follows from (jl2p and the Theorem [T]J^i) . Finally, from the 
proof of the Theorem [T]J^i) it follows that 



E 



„a/(e,)-CG(e,). 



J {eg 



E 



,a/{e,)-CG(e,). 



;J(G(eg)) A;ip(j(G(eg)))^A^, 



which completes the proof of psp in a view of (|12p . 
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4.2 Proof of Theorem [4] 

For a general matrix A with the distinct eigenvalues Aj (hence with the 
independent eigenvectors Sj), under assumption that ql + A is invertible, 
using Prullani integral and the representation A = S diag{Aj} with 
S = (si, . . . , s^r), we can derive the following identity: 

q {ql + A)-^ = exp " ■ (^7) 

Lemma 7. Under assumption |j7p for ^ > 0, 

= expJ / / {exp{-Ct + iax} -l)-e~'^^ P{X{t) £ dx;J{t)) dt\ 

[Jo J [0,00) t J 

•expJ / / {eyip{-Ct + iax} -l)-e-'>^ P{X{t) edx]J{t)) dt\ . 

[Jo J (-00,0) t J 

Proof. By additivity of process X{t) there exists a matrix F such that 
Eiexp{iaX{t)} = exp{F(ia)f} (see Prop. XI. 2. 2, p. 311 of Asmussen 
(2003)). Note that this matrix also has distinct eigenvalues. Then we have: 

/>oo 

J. giaX(e,)-5e,.j(g^) =/ ge"^* exp {- (^1 - F(ia) )t} dt 
L -I Jo 

= exp 1^ (exp {-(a - F{ia))t} - I) ^e"'* dtj 
= exp|y" y (exp{-^t + iaj;} - 1) ^e~^* P(X(t) G dx; J(t)) (itj . 
Note that by identity (I17p the matrices 

/ / {ex.p{-^t + iax} -l)-e-''^ P{X{t) £ dx;J{t)) dt 

Jo J\0,oo) t 



and 



[ [ {exp{-Ct + iax} -l)-e-'^^ P{X{t) £ dx;J{t)) dt 

Jo J(-oo,0) * 



'(-00,0) 

commutes. This gives the assertion of the lemma by the factorization. □ 
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From Lemma [7] and Theorem [T] we have 



exp 



\[ [ {exp{-^t + iax} -l)-e-''^P{X{t) £dx;J{t)) dt\ 
[Jo J[0,oo) t J 

•expJ / / {exY>{-it + iax] -l)-e~'^^V{X{t) edx]J{t)) dt\ 

[Jo J (-oc.fi) t J 

= H(a,OT(a,0, (28) 
where 

II . .. e) = E re^"5(e.)-«G(e,). (29) 



and 



T 

A, 



T(a, = H(g, 0)-Miag (H(g, 0)6^) A'^E e*"^^^')-^^^''^); J(G(eg)) 

(30) 

From Theorem [U^i) it follows that matrices H(a,^) and T(q!,^) are invert- 
ible. Thus, 

H~"'^(a, ^) exp < / / (exp {— + iax} — 1) 

[^0 ^[0,oo) 

^e^*?* P(X(t) G dx; J{t)) dtj 
= T(q, ^) exp i - / / (exp{-^i + zax} - 1) 

[ Jo J{-oofl) 

^6-9* P(X(t) G dx; J(t)) dtj . (31) 

Moreover, each entry of the matrix H(a, ^) is analytical in the upper half 
of the complex plane. The same concerns also the matrix H~^(a,^). Thus 
each entry of the LHS of ()3ip extends analytically to the lower half of the 
complex plane in a and similarly each entry of the matrix on the RHS 
of ()3ip extends analytically to the upper half of the complex plane in a. 
Hence matrices on both sides of (I3ip can be defined in the whole a-plane. 
Observe that each entry of these matrices is bounded function. Indeed, from 
definitions (j29l) and (j30p by Jensen inequality it follows that each entry of 
matrices H(a,^) and T(a,^) is bounded in respective regions. Note that 
reciprocal of determinant of H(a,^) is also bounded within any bounded 
circle. Thus on any circle each entry of H~^(a,^) is bounded. Similarly, 
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one can prove that each entry of the second factors of the RHS and LHS of 
(|3ip is bounded. Thus by Liouvihe's Theorem each entry of (j3ip must be a 
constant. Putting a = = gives the assertion of the theorem. 



4.3 Proof of Theorem m 

We will use now the martingale technique to prove Theorem [5] (see Borovkov 
and Burq (2001) for similar considerations in the case of spectrally negative 
Levy processes). 

Lemma 8. We have, 

/ ^^^^J^ = / te-"* / — y ^45Mp,(r+ G dt; Jit) = j). 

Jy ^'{Hq)) Jo Jv ^ MHq)) 



y i^'mq)) Jo 

Proof. From ([4]) we have 

e-*^- =E,e-''-^^h^(^+)($(g))/hi($(g)) . 
Differentiating with respect to q and noting that ^'{q) = 1/ n' {^(q)) we have 

n'i^iq)) Jo M^iQ)) 

The proof completes by dividing left-hand side and right-hand side of ([32|) 
by X and integrating them w.r.t. dx over (y, oo). □ 



Lemma 9. We have, 

N 



r^dxe-^^ .oo /-oo " h,(c|>(g)) 

Jy ^'m)) Jo Jy friMm) 



{X{t) G dx; J{t) = j) dt. 



a = «'($((z))Ef + h,($(9)) - Ef^''\.^.mq)) 



Proof. From Corollary XI. 2. 6 of Asmussen (2003) we have 

'j(r+) 

and hence 

E*('?V+ 1 
lim -1 ^ = . (33) 

a->(» a K'[^{q)) 

Let Ta = J^lA{y^it)) dt be the time spend by our process in the set A. 
Note that ^ 

T+ - T(„^,o] < T(o,a] < T+ + T^l^^'^Jj , (34) 
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where T^l^^j denotes the time spend in (— cxd,0] by the process {{X{t + 
r+) - a, J(r+ + t)),t > 0}. Moreover, 

maxEf(^)r(_o,,o] = maxE* ("^^ / I(_^,o](^(t)) dt 

« ' ' ^ Jo 

POO 

» Jo 

< maxEf" / e-""'> dt < max / Ef"'e-»^"> dt 



id hi{^{q)) Jo q id h,($(g)) 

Taking expectation from both sides of ()34p and dividing by a we derive: 

hm = . (35) 

Observe that for 0<a<6<c<oo 

and hence 



which together with (j35p gives: 



^'(4.(9)) 

Above can be rewritten in the following way: 

f«)(X(t)G(0,a]) 

N 







J{0,a] MHl)) 



J{0M Jo 
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Since this equation holds for any a > 0, we have 

N 



dx 



-qt 



hj(^(g)) 



fri MHq)) 



Fi {X{t) G dx; J{t) = j) dt 



which completes the proof by integrating them w.r.t. dx over (y, oo). □ 
From Lemma [8] and [9] we have the following equality of Laplace trans- 
forms: 



te 



-qt 



dx sr^ \ij{<^{q)) 

y X Mm) 



(r+ G dt; J{t) = j) 



-qt 







N 



that implies the equality of the following measures: 

N 



^h,(cl>(g)) (t— P.(t+ G dt; J{t) = j) - Fi (Xit) G dx; J{t) = j) dt) = 
7 = 1 ^ ^ ^ 



for each i = \^ . . . ,N . This is equivalent to: 

(tdxP(r+ G dt; J{t)) - xP G dx; J{t)) dt) h($(g)) = . 



(37) 



Choosing now the independent vectors h(<I>(gi)), h(<^(g2)), • • • , h($(gAr)) com- 
pletes the proof. 

4.4 Proof of Theorem [6] 

By Kendall's identity given in Theorem [5] it suffices to prove that 

P(X(t) G dx, I{t) = 0; J{t))dt = ^^ij^ e dt; J{t))dx 
or that for all (7 > and sufficiently large s > 0: 

/•oo /'OO 

g / e-9*dt / e^^P(X(t) G (ix,/(t) = 0; J(t)) 



c Jo 



/"OO /"OO 

/ e-«* / e"^(ixP(r+ G dt; J(t)) 
JO JO 



(38) 



that is equivalent to 



hm E 

a— +00 



^sX(e,)+a/(e,). 



]J{^q 



e-«^- ; J(t+) . (39) 
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We prove (j39p passing from its left-hand side to its right-hand side. Let 
q = q — k{s). The change of measure (jl]) and Wiener-Hopf factorization 
given in Theorem [TJii) yields: 



lim E 

a— >oo 



lim Ah(s)E'^ 



J(e,) Ah(s)-i 



lim ^Ah(s)E^ e°^(^);J(e^) A^{s)-' 
Q— >oo g L 



hm Ah(s)g ( - F,(a)/a ) A;3iH,(g, -a)^ H,(g, 0)"^ ^ A,>h(s)''. 



a 



Note that 



and from ([7]) 



lim (^I-F,(a)/a = (-d)"^ 



1 

-I 

c 



lim -E,s{q,-a)'^ 

a^oo a 



I. 



Hence 
hm E 

Q— >00 



-Ah(s)A: 



H,(g,0)-i A^^Ah(s)-^ 



Using classical arguments for reversed process we can proceed as follows: 



c L J 



-Au(s) 



e~~^^^-J{Tt) A^^dxA^is)-' 



AJs)B' 



J{t+) Ah(s)~^ dx 



c Jo 



E 



-(5+K(s))rj' +SX. J^^+ 



dx 



which is the right-hand side of ()39p . 



Acknowledgements 

ZP would like to thank Andreas Kyprianou who works with him at the 
beginning of this project and who gave lots of valuable comments. This work 
is partially supported by the Ministry of Science and Higher Education of 
Poland under the grant N N2014079 33 (2007-2009). 



17 



References 

[1] Asmussen, S. (2003) Applied probability and queues. Sec. Edition. 
Springer- Verlag, New York. 

[2] Asmussen, S. and Kella, O. (2000) A multi-dimensional martingale for 
Markov additive processes and its applications. Adv. in Appl. Probab. 
32(2), 376-393. 

[3] Arjas, E. and Speed, T.P. (1973) Symmetric Wiener-Hopf factoriza- 
tion in Markov additive processes. Z. Wahrscheinlichtkeittheoriee verw. 
Gebiete 26, 105-118. 

[4] Bertoin, J. (1996) Levy processes. Cambridge University Press. 

[5] Bertoin, J. (2000) Subordinators, Levy processes with no negative 
jumps, and branching processes. Notes from the course given at Ma- 
PhySto, Denmark. 

[6] Borovkov, K. and Burq, Z. (2001) Kendall's identity for the first cross- 
ing time revisited. Electronic Communications in Probability 6, 91-94. 

[7] Breuer, L. (2008) First passage times for Markov-additive processes 
with positive jumps of phase-type. J. Appl. Prob. 45(3), 779-799. 

[8] Dicker, T. and Mandjes, M. (2006) Extremes of Markov-additive pro- 
cesses with one-sided jumps, with queueing applications. Preprint 
( )http: / / www.proba.ucc. ie/^td3 / publications / modulatedfjuid.pdf, ) . 

[9] Feller, W. (1970) An Introduction to Probability Theory and Its Appli- 
cations, Vol. II. 2nd edition, John Wiley, New York. 

[10] Fristedt, B.E. (1974) Sample functions of stochastic processes with sta- 
tionary independent increments. Advances in probability and related 
topics 3, 241—396. Dekker, New York. 

[11] Greenwood, P. and Pitman, J.W. (1980) Fluctuation identities for Levy 
processes and splitting at the maximum. Adv. Appl. Prob. 12, 893-902. 

[12] Ivanovs, J., Boxma, O. and Mandjes, M. (2008) Singularities of the 
generator of a Markov additive process with one-sided jumps. Eurandom 
Report 2008-37. 



18 



[13] Kaspi, H. (1982) On the symmetric Wiener-Hopf factorization for 
Markov additive processes. Z. Wahrscheinlichtkeittheorie verw. Gebi- 
ete 59, 179-196. 

[14] Kyprianou, A. (2006) Introductory Lectures on Fluctuations of Levy 
Processes with Applications. Universitext, Springer. 

[15] Kyprianou, A.E. and Palmowski, Z. (2008) Fluctuations of spectrally 
negative Markov additive processes. Sem. de Probab. XLI, 121-135. 

[16] Maissoneuve, B. (1975) Exit systems. Annals of Probability 3, 399-411. 

[17] Pistorius, M.R. (2005) A potential theoretical review of some exit prob- 
lems of spectrally negative Levy processes. Sem. de Probab. 38, pp. 
30-41. 

[18] Prabhu, N. U. (1998) Stochastic storage processes. Queues, insurance 
risk, dams, and data communication. Sec. Edition. Springer- Verlag, 
New York. 

[19] Presman, E.L. (1969) Factorization methods and boundary problems 
for sums of random variables given on Markov chains. Math. Izv. Akad. 
Nauk USSR 3, 815-852. 

[20] Rogers, L.C.G. (1994) Fluid models in queueing theory and Wiener- 
Hopf factorization of Markov chains. Ann. Appl. Probab. 4, 390-413. 

[21] Rogozin, B.A. (1966) On the distribution of functionals related to 
boundary problems for processes with independent increments. Theory 
Probab. Appl. 11, 580-591. 

[22] Spitzer, F. (1964) Principles of Random Walk. New York, Van Nos- 
trand. 



19 



